Multiphoton ghost imaging with classical light 
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One of the possible types of n-th order ghost imaging is experimentally performed using multi-photon 
(higher-order) intensity correlations of pseudothermal light. It is shown that although increasing the 
order of intensity correlations leads to the growth of ghost imaging visibility, it at the same time reduces 
the signal-to- noise ratio. Therefore, ghost imaging with thermal light is optimal in the second order in 
the intensity. 
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INTRODUCTION 

The technique of ghost imaging (GI), first demon- 
strated in 1995 [1 j, still attracts much interest in connec- 
tion with remote sensing [2 , lenseless imaging [3 , and 
fundamental problems, such as the boundary between 
quantum and classical physics (see, for instance, [4 j). The 
idea of ghost imaging is based on the correlations exist- 
ing between the intensities of two beams, which will be 
further called the signal one and the reference one. The 
object to be imaged (a mask) is placed into the signal 
beam, and the photons transmitted through it are reg- 
istered by a "bucket" detector, which collects all signal 
radiation and hence does not resolve the shape of the 
mask. The detector in the reference beam, on the con- 
trary, is spatially resolving, but there is no mask in the 
reference beam. The coincidences, or photocurrent cor- 
relations, of signal and reference detectors are registered 
as a function of the reference detector position. When 
the reference detector registers photons whose correlated 
signal counterparts pass through the mask, there is an 
increase in the coincidence counting rate (photocurrent 
correlation function). Thus, by scanning the reference 
detector one restores the image of the mask in the signal 
beam. 

Although ghost imaging was at first supposed to re- 
quire two-photon light generated via spontaneous para- 
metric down-conversion (SPDC), it was later demon- 
strated with classical light [T7] , in particular, pseudo- 
thermal light obtained by passing a coherent beam 
through a rotating ground-glass disk |6j [7J [8] and true 
thermal light [9]. In a large number of works, both near- 
field and far-field ghost imaging with pseudo-thermal 
light was performed, as well as ghost imaging of a purely 
phase object [10 or a reflecting object [2]. It seemed 
that the only advantage of using nonclassical light for 
ghost imaging was a higher visibility than in the case of 
thermal light. However, in several works it was proposed 
to increase the visibility of ghost imaging with classical 
light by increasing the value of normalized correlation 
function. In Ref. [11, it was performed by introducing a 
threshold for the bucket detector output and hence choos- 



ing only peaks of the intensity fluctuations. In Ref. [12], 
it was suggested to use the fact that the normalized in- 
tensity correlation function (CF) of the n-th order for 
thermal light grows as n\ and thus to increase the visi- 
bility of ghost imaging with thermal light by passing to 
higher-order intensity correlations. This idea was real- 
ized recently [13 and the growth of ghost-imaging visi- 
bility with the order of the registered intensity moment 
was indeed demonstrated in experiment. 

However, as it was pointed out in Refs. [12] [14] [17] , it is 
not the visibility of a ghost image that matters but rather 
its signal-to-noise ratio (SNR) . The noise of a ghost image 
is related to the variance of the intensity moment under 
measurement and contains both the shot-noise compo- 
nent and the excess noise component, as well as inter- 
mediate terms [14]. Theoretical quantum calculation of 
the noise for ghost imaging with thermal and nonclassi- 
cal light has been performed in Ref. [13] ; in Ref. [15] the 
noise for ghost imaging with pseudo-thermal light was 
also measured experimentally. In Ref. [17 the so called 
"contrast-to-noise ratio" has been theoretically studied 
for high-order thermal ghost imaging using a classical 
approach. However, the present work is the first consid- 
eration and experimental measurement of SNR for mul- 
tiphoton ghost imaging with pseudo-thermal light, i.e., 
ghost imaging based on higher-order intensity moments. 



THEORY 

There are many ways to generalize ghost imaging in 
terms of using higher order intensity correlations. First 
for the sake of simplicity we consider an (n-l)-port 
Hanbury-Brown and Twiss (HBT) setup depicted in Fig. 
1. As in a standard lensless ghost imaging setup, radia- 
tion from a thermal source (usually a rotating ground- 
glass disk illuminated by a laser beam) is split on a 
beamsplitter (BSi). One of the two output beams (the 
reference one) is further split on n-2 beamsplitters and 
detected in the far- field zone by n-1 spatially resolving 
detectors at locations xi,...,x n _i. The other (signal) 
beam is sent to a bucket detector (B) with a mask (M) 
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beam transmitted through the mask, 1 < k < M . the 
correlation function [3] will take its maximum value, 



Gtl = {n-l)\{M + n-l)I n . 



(5) 
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Otherwise, if k is outside the interval [1,M], the cor- 
relation function will take the 'background' value, 



FIG. 1: Schematic of an n-port HBT setup for measuring an 
n-point n-th order CF. 



placed in front of it. Using all n detectors, one measures 
the intensity correlation function of order n 
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The visibility of nth-order ghost imaging can be de- 
fined as 
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where I(x^) and J(y) are instant intensities at corre- 
sponding locations measured by (n-1) reference and one 
signal detectors, respectively. Angle brackets correspond 
to time-averaging and transmission of the mask is de- 
noted by the function T(y) taking only values or 1. 
Because the bucket detector collects all signal radiation, 
integration is implied over its area A s . Since the bucket 
detector is supposed to be larger than the mask, integra- 
tion can be done over the mask area. 

Since it is rather complicated to analyze properties of 
an n-point correlation function (or (n-l)-point ghost im- 
age) we focus on a more simple 2-point case where posi- 
tions of all (n-1) reference detectors coincide: x\ = . . . = 
x n -\ — x. Eq. ([I]) then turns into 

G (n) (x) ^ (/ n-l (x) f J (y)T(y)dy) . (2) 

J A s 

To simplify the calculation, we will consider one- 
dimensional case and use a very convenient approach of 
Ref. [15] where a discrete set of modes was introduced. It 
means that from continuous intensity distributions /(x) 
we pass to the values of intensities at discrete points X&, 
and integration in Eq. Q is replaced by summation over 
the modes k' transmitted by the mask, 1 < k' < M: 
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Calculation of the intensity correlators from Eq. (|3| 
can be easily performed taking into account that for ther- 
mal light 

(AW) = {( l + ™) ! Skk' + U ml (1 - 6 kk ,)}I l+m . (4) 

Here, / is the time-averaged intensity, which is assumed 
to be uniform over the whole area of the mask, and 5k,k' 
is the Kronecker symbol. Hence, for the positions of the 
reference detector that have counterparts in the signal 
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i.e., the visibility grows with the order of the correlation 
function measured. From Eq.(j8|, one can see that the 
visibility decrease due to a large number M of modes in 
the image can be reduced by passing to a high order n. In 
the single- mode case, M = 1, the values of visibility for 
different orders are = 1/3, = 1/2, V^ 4 ) = 3/5,... 
Note that the results for the maximum and background 
values of the correlation function, as well as for the visi- 
bility, will be the same if quantum-mechanical approach 
is used in the calculation. However, in the calculation of 
the noise it is important to use the quantum-mechanical 
treatment (as, for instance, in Ref. [14], since it is the only 
way to take into account the shot noise. For this purpose, 
we will write the expression for the maximum and back- 
ground correlation functions in the quantum form, 
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Here, angular brackets denote averaging over the thermal 
state of M + 1 modes whose photon creation and anni- 
hilation operators are a^, a^, k = 0, 1, . . . M. The mode 
k = is the only one considered outside of the mask. 
We will now calculate the signal-to-noise ratio, taking 

as the signal the difference S = Gmix — G^aLk anc ^ as tne 
noise, the standard deviation of this value, which is 



AS = V Var(G^2 x ) + Var(G^ ck ) - 2Cov(G : 



(n) 



(n) r {ri) x 
^back/' 

(io) 

where Var(Gm2x) and Var(G^ ck ) are variances of the 
maximum and background correlation functions and 
Cov(G^ x ,G[ ) " ) ck ) is their covariance. The latter has to 
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be taken into account since Gmax,G|[>ack are not statisti- 
cally independent. 



Calculation of the values entering Eq.(lO) is done by 



passing to normally-ordered correlators using Wick's the- 
orem [16]. For instance, Var(G^ ck ) is calculated as 

Var(G^ k ) = 

M M 

= E E ({alY-'ala k ar\alY-\l,a k ,ar 1 ) 

k=l k' = l 

-[(n-l)!] 2 M 2 / 2n , (11) 



which yields 
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Similarly, calculation of Var(Gmix) 
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and 



Cov(Gmax, G^ k ) leads to the following expression 
for the signal-to-noise ratio: 
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In the low-intensity limit, this expression becomes 
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We see that the signal-to-noise ratio grows polynomi- 
ally with the intensity of light, the growth index being 
n/2. By plotting dependencies ( p~3| ) for various orders n 
for single-mode case (Fig. 2), we see that although the 
growth is faster for higher n, the signal-to- noise ratio sat- 
urates at large intensities, and the maximal values of SNR 
at higher orders are smaller. Indeed, the high-intensity 
asymptotic values of SNR are 
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i.e., they decrease at large n. Analogous calculation 
for the case of SPDC gives the following expression for 
SNR 
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which coincides with the result obtained for second-order 
thermal ghost imaging. 
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where m is the mean number of photons per mode. In 
a high-intensity limit 



FIG. 2: Normalized SNR dependence on intensity expressed 
in photons per mode for 2-nd (red solid), 3-rd (blue dashed) 
and 4-th (orange dot-dashed) order ghost images for single- 
mode case (M — 1). 



Fig. 3 shows the dependence of the normalized SNR 
on the intensity expressed in photons per mode for SPDC 
and second-order thermal light ghost imaging for single- 
mode and M = 10 cases. As one can see SPDC results 
in higher SNR values only at low intensities, reaching 
a maximum value of about 0.27 at intensity of approx- 
imately 0.8 photons per mode for the single-mode case. 
This maximum is just 6% higher that the hight-intensity 
limit (which is ^7= w 0.26), so for single- mode case 
the difference between the maximum SNR obtained with 
biphotons and the one obtained with thermal light is neg- 
ligible. For M = 10 SPDC has a maximum SNR value 
of 0.11 at 0.07 photons per mode which is 90% higher 



4 



0,1 



0,01 



1E-3 



1E-3 0,01 0,1 



FIG. 3: Normalized SNR dependence on intensity expressed 
in photons per mode for single-mode SPDC ghost imaging 
(green solid), thermal light GI (red dashed), SPDC ghost 
imaging for M = 10 (purple dash-dotted) and thermal light 
GI for M = 10 (blue dotted). Given in the inset is the central 
part plotted in linear scale. 



that the high-intensity limit (which is 



0.06 for 



this case). 

Experimentally it is more convenient to measure a 2- 
point correlation function using a 2-port HBT setup (Fig. 
4) instead of an n-port setup. Below, we provide a proof 
of equivalence of the above mentioned setups for strong 
thermal light (which is the case of consideration) in terms 
of measuring a 2-point n-th order intensity CF. In quan- 
tum description, an n-point n-th order CF is denned by 
the following expression 
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where a\ and hi are the photon creation and annihilation 
operators at point Xi. If n-1 spatial arguments of the CF 
coincide x\ = x<z = . . . = x n _i, it becomes 
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The coincidence counting rate measured by the setup 
shown in Fig. 4 corresponds to a different 2-point CF, 
having the form 



G (n) (xi,xi,...,xi,x n ) = la\a 1 a\a 1 ...a\ l a n 



(20) 



Using the commutation relations one can represent CF 



( 20 ) as a series of normally-ordered CFs [T6] 
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2 = 2 

with C n = 1. Each CF C?W contains 2i-2 operators taken 
at point x\ and two operators taken at point x n . 

Since for thermal light (?W scales as the i-th order of 
intensity, in the high-intensity limit the term with i = n 



in (21) dominates, and CFs (19) and (20) coincide. In 



this case, our measurement procedure gives the same as 
the n-port HBT setup. 
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FIG. 4: Schematic of a 2-port HBT setup for measuring a 
2-point n-th order ghost image. 
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FIG. 5: Experimental setup for measuring n-th order ghost 
images. 



EXPERIMENT 

Experimental setup is presented in Fig. 5. As the ra- 
diation source, we use a frequency doubled Q-switched 
Nd:YAG laser with the wavelength 532 nm, pulse dura- 
tion 5 ns, and the repetition rate of 47 Hz. The laser 
beam with diameter set by aperture Ai is projected onto 
a rotating ground-glass disk GGD to form pseudothermal 
light. The central part of the speckle field is selected by 
aperture A 2 . Polarizing filter P is introduced to select 
a linear polarization, thus eliminating any polarization 
effects introduced by the mirror. Reflecting from mirrors 
Ml and M2 the beam is divided into two by a nonpo- 
larizing beam-splitter BS. Transmitted (signal) and re- 
flected (reference) beams are reflected by mirrors M3 
and M4 respectively. 30-degrees prism PR makes both 
beams parallel and sends them to a digital photo camera 
DPC. Lens L is used to image the desired area of the 
speckle field onto the photo camera's matrix 20.7x13.8 
mm (2640x1760 px). 

In experiment it is more convenient to use two differ- 
ent parts (signal and reference) of a single CCD camera 
to capture both signal and reference beams instead of us- 
ing a CCD and a bucket detector as there is no need to 
synchronize the latter two. In this case the signal from 
the bucket detector is obtained by calculating the total 
signal captured by the reference part of the CCD cam- 
era. The idea of using only one CCD camera has been 
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FIG. 6: Single captured frame. The area marked by the solid 
and dashed white rectangles represent two correlated parts of 
the speckle field that are selected for ghost image calculation. 

recently implemented in experiment [15] . In our case the 
"software" mask was calculated using a true "physical" 
realization of the reference beam instead of using only 
one beam as both signal and reference. 

Exposure time was chosen to be 1/50 seconds, so that 
only one pulse was captured at every frame and the frame 
rate was about 0.2 fps. Ground glass disk rotation speed 
was small enough to consider it stationary for the dura- 
tion of laser pulse and large enough to allow considerable 
(compared to GGD surface inhomogeneity size) linear 
displacement between two consequent frames. Thus all 
the captured frames were time independent and a total 
of 5000 frames were captured. A typical frame is pre- 
sented in Fig. 6. The average speckle size was about 30 
px (pixels) which correspond to 150/i. By summing the 
intensities over a desired area (i.e. mask) we obtained the 
"bucket detector" signal which was correlated with every 
pixel's signal in the signal arm. Two correlated parts of 
the signal and reference arms were used for ghost image 
calculation - they are depicted by solid and dashed white 
rectangles in Fig. 6. 

For clarity we choose a very simple object to image 
- a single slit represented by a line of 1 px height and 
variable width. By changing the length of the line one can 
change the number of modes M. An array of signal values 
S(i, j, k, I) is calculated as the difference of Gmi x (i, j) and 
^back(^> 0- Mean value S of j, Z) is then calculated 
as follows 

S = (GgL(i,j)) M - (G%L(k,l)) B , 

where angle brackets denote averaging that is done over 
areas Ai and B. The first area is a line along the the 
slit with j = 2 slit (vertical position of the slit), the lat- 
ter is a rectangular area of a ghost image outside of the 
slit. Noise is calculated as the standard deviation of 
S(iJ sHt ) = (SM(iJ sHu k 1 l)) B . 

As one can see from Fig. 7, SNR (the "clarity" of 
the slit's ghost image) decreases with the growth of the 
ghost image order as it is predicted by the theory. For 
a bigger object (a longer slit) this decrease is even more 
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FIG. 7: Ghost images of a short and long one-dimensional 
single slits for 2-nd (a), 3-rd (b) and 4-th (c) order ghost 
imaging measured using 5000 frames. 
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FIG. 8: Visibility dependence on parameter M (number 
of modes) for the 2-nd (red), the 3-rd (blue) and the 4-th 
(orange) order ghost imaging of a single slit. 

pronounced: the longer slit image has almost drowned in 
the noise for the case of fourth-order ghost image. By 
changing the size of the slit in steps one can plot the 
dependence of visibility and SNR on the size of the mask. 

From Fig. 8 it is seen that visibility increases for 
higher-order ghost images and the values obtained are 
in a good agreement with the theoretical fit (solid lines) . 
Due to laser power instability over time and GGD inho- 
mogeneity, the fit is made by assuming that the observed 
values of autocorrelation functions differ from the theo- 
retical values by constant factors = F n n\ with their 
values close to unity. As far as the "real" speckle size 
might be different from the one determined as a FWHM, 
the parameter M is assumed to be different from a the- 
oretical one by a fixed value S. F n were independently 
calculated as mean values of autocorrelation functions of 
the corresponding order over 5 million pixel areas and 
make up the following values for second-fourth orders: 
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FIG. 9: Normalized SNR dependence on the parameter M 
(number of modes) for the 2-nd (red), the 3-rd (blue) and the 
4-th (orange)order ghost imaging of a single slit. 



F 2 = 0.988, F 3 = 0.984 and F A = 0.997. The best fit 
corresponded to S = 1.19, which means that the num- 
ber of modes is determined by speckle width (FWHM) 
multiplied by S. 

In experiment one can only measure sample mean and 
sample standard deviation, so, in fact, all SNR values are 
multiplied by a factor of y/N, where N is the total number 
of images used to calculate the ghost image. Fig. 9 shows 
the dependence of the normalized SNR on the number of 
modes M. One can see that although the visibility does 
increase with the order, the SNR drops. In order to have 
the same SNR together with higher visibility values one 
should increase the number of images TV used. 



CONCLUSION 

We have experimentally performed the 2-nd, the 3-rd 
and the 4-th order 2-point ghost imaging in pseudother- 
mal light using a 2-port HBT setup and theoretically 
proved that such measurements are equivalent to those 
made with 2-, 3- and 4-port setups respectively. We 
showed that for higher order ghost images defined as in 
expression ([2| there is indeed an increase in imaging vis- 
ibility and decrease in SNR at the same time. It is found 
that the second order ghost imaging is optimal for ther- 
mal light in the sense of the highest SNR values compared 
to those in higher orders for a fixed sample size (number 
of images). Higher visibility values might be obtained us- 
ing higher order correlations but at the expense of SNR 
decrease which, in turn, can be compensated by increas- 
ing the sampling size. 

In the single-mode case thermal light has slightly lower 
SNR values than SPDC (the difference is less than 6%). 
In multi-mode case, the difference is larger: ghost imag- 
ing with SPDC has a higher SNR maximum compared to 
thermal light. This result is in agreement with theoretical 
consideration in [14] . But at high intensities SNR values 



are the same for both cases at any values of M. From 
this point of view, if one is interested in best SNR values, 
thermal light and SPDC give identical result, whereas if 
one is interested in having the best visibility at a fixed 
number of images used SPDC has an advantage. 

Further research is needed to investigate other cases 
of ghost imaging experimentally. For instance, recent 
results p~3j [17] suggest that higher SNR can be obtained 
by exploiting n-th order ghost imaging with symmetrical 
orders of reference and signal intensities. 
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